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(54) Ultrasound contrast imaging with double-pulse excitation waveforms 



(57) The method of imaging the nonlinear compo- 
nents of ultrasound echo signals returned from scatter- 
ed comprises the steps of constructing a double-pulse 
excitation waveform, the two pulses composing this 
waveform having spectra with known relative frequency 
dependencies, differing in amplitude and phase so that 
the double-pulse waveform can be considered as a con- 
volution of any single of said pulses with a known coding 
function, transmitting said double-pulse waveform with- 
in a single transmission of an ultrasound beam in the 



direction of said scatterers, deconvolving the received 
rf-waveform by a first appropriate decoding function to 
obtain a new rf-waveform, deconvolving said received 
rf-waveform by a second and different appropriate de- 
coding function to obtain another rf-waveform, realign- 
ing in time both deconvolved rf- waveforms obtained 
from the same returned echo signal by a time delay as 
determined by the transmit-coding function, normalizing 
amplitudes and taking the sum of both rf-waveforms to 
effectively cancel all echo components caused by linear 
scatterers. 
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Description 

[0001] The invention relates to a method of ultrasound 
imaging of organs and tissues by detection of ultrasound 
backscatter from a region which may contain nonlinear 
scatterers such as microbubbles used as a contrast 
agent, the method comprising projecting an ultrasound 
beam to a zone of tissue to be imaged, receiving the 
echo reflected from the tissue as a radiofrequency re- 
sponse signal, processing the radiofrequency response 
into a demodulated video output signal, storing the out- 
put in a video scan converter, and scanning the tissue 
to produce a video image of the region under investiga- 
tion. The invention also comprises a system for ultra- 
sonic imaging of organs or tissues which may contain 
nonlinear scatterers such as microbubbles used as a 
contrast agent, the system comprising an ultrasonic 
probe for transmitting and receiving ultrasonic signals, 
signal processing means, means for storing the proc- 
essed signals and a display element. Use of the system 
for imaging of organs, tissues and blood vessels is also 
disclosed. 

Background Art 

[0002] Wide acceptance of ultrasound as an inexpen- 
sive non-invasive diagnostic technique coupled with 
rapid development of electronics and related technology 
has brought about numerous improvements to ultra- 
sound equipment and ultrasound signal processing cir- 
cuitry. Ultrasound scanners designed for medical or oth- 
er uses have become cheaper, easier to use, more com- 
pact, more sophisticated and more powerful instru- 
ments. However, the changes of acoustic impedance 
occurring within the living tissue are small and the ab- 
sorption of ultrasound energy by different types of tissue 
(blood vessels, organs, etc.) are such that some diag- 
nostic applications remain unmet challenges, despite 
these technical developments. This situation changed 
considerably with the development and introduction of 
administrate ultrasound contrast agents. Introduction 
of contrast agents made from stabilized suspensions of 
gas microbubbles or microballoons into the bloodstream 
and organs to be investigated have demonstrated that 
better and more useful ultrasound images of organs and 
surrounding tissue may be obtained with ultrasound 
equipment. Thus, pathologies In organs like the liver, 
spleen, kidneys, heart or other soft tissues are becom- 
ing more readily recognizable, opening up new diagnos- 
tic areas for both B-mode and Doppler ultrasound and 
broadening the use of ultrasound as a diagnostic tool. 
[0003] Recently, ultrasound techniques, i.e. scan- 
ners, electronic circuitry, transducers and other hard- 
ware and software components are showing great 
progresses in their abilities to exploit, to a fuller extent, 
the specific properties of ultrasound contrast agents. 
This is made possible by the vision by ultrasound instru- 
ment manufacturers of the vast potential offered by 



these contrast agents towards more accurate diagnosis, 
thanks to enhanced imaging capabilities and quantifica- 
tion of blood flow and perfusion .Thus, what was almost 
independent developments of these related segments 

5 of the field are now providing opportunity to draw on syn- 
ergies offered by studies in which the electronic/ ultra- 
sound characteristics of the apparatus and the physical 
properties of the contrast agent are combined. A few ex- 
amples of such studies reported improvements from 

10 specific agents/equipment combinations, such as har- 
monic contrast imaging. These synergies are thus open- 
ing new areas of experimentation, innovation, and 
search of more universal methods for producing greater 
tissue resolution, better image and greater versatility of 

is ultrasound as a diagnostic technique. There is no doubt 
that, provided their implementation is kept relatively sim- 
ple, these will be widely accepted. 
[0004] An attempt towards improved ultrasound im- 
aging is described in WO-A-93/1 2720 (Monaghan) 

20 which discloses a method of imaging of a region of the 
body based on subtracting ultrasound images obtained 
prior to injection of a contrast agent from the images of 
the same region obtained following administration of the 
contrast agent. Based on this response subtraction prin- 
ts ciple, the method performs superposition of images ob- 
tained from the same region prior to and after adminis- 
tration of the contrast agent, providing an image of the 
region perfused by the contrast agent freed from back- 
ground image, noise or parasites. In theory, the method 

30 described is capable of providing a good quality images 
with enhanced contrast. However, in practice, it requires 
maintenance of the same reference position of the re- 
gion imaged for a long period of time, i.e. long enough 
to allow injection and perfusion of the contrast agent and 

35 maintenance of an enormous amount of data. Therefore 
the practical implementation of the method is very diffi- 
cult if not impossible. The difficulty is partly due to inev- 
itable internal body movements related to breathing, di- 
gestion and heart beat, and partly due to movements of 

40 the imaging probe by the ultrasound operator. Most re- 
altime imaging probes are commonly handheld for best 
perception, feedback and diagnosis. 
[0005] Interesting proposals for improved imaging of 
tissue containing microbubble suspensions as contrast 

45 agent have been made by Bums, P., Radiology 185 P 
(1992) 142 and Schrope, B. et al., Ultrasound in Med. 
& Biol. 19 (1993) 567. There, it is suggested that the 
second harmonic frequencies generated by non-linear 
oscillation of microbubbles be used as imaging param- 

so eters. The method proposed is based on the fact that 
normal tissue does not display nonlinear responses to 
the same extent as microbubbles, and therefore the sec- 
ond harmonics method allows for contrast enhance- 
ment between the tissues with and without contrast 

55 agent. Although attractive, the method has its shortcom- 
ings, as its application imposes several strict require- 
ments. Firstly, excitation of the fundamental "bubble- 
resonance" frequency must be achieved by fairly nar- 
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row-band pulses, i.e. relatively long tone bursts of sev- 
eral radio-frequency cycles. While this requirement is 
compatible with the circuits and conditions required by 
Doppler processing, it becomes hardly applicable in the 
case of B-mode imaging, where the ultrasound pulses 
must be of very short duration, typically one-half or one- 
cycle excitation. In this case, insufficient energy is con- 
verted from the fundamental frequency to its "second- 
harmonic", and thus the B-mode imaging mode cannot 
greatly benefit from this echo-enhancing method. Sec- 
ondly, the second harmonic generated is attenuated, as 
the ultrasound echo propagates in tissue on its way back 
to the transducer, at a rate as determined by its frequen- 
cy, i.e. at a rate significantly higher than the attenuation 
rate of the fundamental frequency. This constraint is a 
drawback of the "harmonic-imaging" method, which is 
thus limited to propagation depths compatible with ul- 
trasound attenuation at the high "second-harmonic" fre- 
quency. Furthermore, in order to generate echo-signal 
components at twice the fundamental frequency, "har- 
monic imaging" requires non-linear oscillation of the 
contrast agent. Such behavior imposes the ultrasound 
excitation level to exceed a certain acoustic threshold 
at the point of imaging (i.e. at a certain depth in tissue). 
During nonlinear oscillation, a frequency conversion 
takes place, causing in particular part of the acoustic 
energy to be converted from the fundamental excitation 
frequency up to its second harmonic. On the other hand, 
that level should not exceed the microbubble burst level 
at which the microbubbles are destroyed, and hence 
harmonic imaging will fail due to the destruction of the 
contrast agent in the imaging volume. The above con- 
straints require that the imaging-instrument be set-up in 
such a way as to ensure the transmit-acoustic level to 
fall within a certain energy band: high enough to gener- 
ate second harmonic components, but low enough to 
avoid microbubble destruction within a few cycles. 
[0006] Thus, methods which treat electronic echo sig- 
nals during normal realtime ("on the fly") investigations 
are those most desirable, allowing better imaging and 
wider use of ultrasound diagnostic imaging. Such meth- 
ods are based on an enhancement of the echoes signals 
received from the regions imaged, using signal process- 
ing functions which are designed to enhance the con- 
trast between regions containing contrast agent from 
those without contrast agent, on the basis of nonlinear 
or frequency-dependent parameters, would be simple 
to use and implement in new instrument designs. 
[0007] US 5'632'277 and US 5706'819 already dis- 
close methods and apparatus for the detection and im- 
aging of harmonic echo components from microbubble- 
based contrast agents in blood. These methods utilize 
first and second ultrasound pulses that are alternatively 
transmitted into the medium being imaged. The first and 
second ultrasound pulses are amplitude modulated sig- 
nals in the radiofrequency range, the first ultrasound 
pulse differing only in sign (polarity) from the second ul- 
trasound pulse transmitted. The echo signals generated 



by these successive pulses are stored in memory and 
combined by adding them so that the linear components 
cancel, leaving only the nonlinear component to be im- 
aged. Accordingly, since tissue generally reflects less 
harmonic components than microbubbles, such 
processing enables the microbubble echoes of the con- 
trast agent to be received with a high signal to noise 
ratio. 

[0008] In these known approaches, for each " line-of- 
sight ", or in otherwordsfor each ultrasound beam steer- 
ing and focusing properties, a minimum of two succes- 
sive pulses is required to cancel echoes from tissues 
while preserving significant signals of echoes from mi- 
crobubbles. The radio-frequency echo signals are 
stored in a temporary memory following excitation of 
each successive pulses. 

[0009] Alternative methods have been proposed in 
US 5'961'463, or WO 99/30617 according to which, by 
alternating the polarity or otherwise coding the succes- 
sive excitation conditions and taking the sum, difference 
or other combinations of these rf-echo signals, essen- 
tially zero signals are produced from linear reflectors 
otherthan the microbubbles, while echo signals from mi- 
crobubbles produce non-linear signals that can be used 
to construct images with great sensitivity and enhanced 
contrast between microbubbles and tissues not contain- 
ing a significant number of microbubbles. 
[0010] There are several disadvantages common to 
these techniques. Firstly, because of the finite propaga- 
tion velocity of sound in tissues, a time interval of several 
hundred microseconds must elapse between succes- 
sive transmit pulses, in order for the echoes from the 
deepest regions of interest to return to the ultrasound 
probe before a following transmit pulse can be applied. 
This requirement limits the achievable frame rate, which 
is reduced in a proportion depending on the number of 
pulses fired in each direction, compared to similar im- 
aging conditions without contrast-specific transmit cod- 
ing. Secondly, the bubbles can move significant distanc- 
es in the blood vessels in the time interval between puls- 
es. Because of that, the resulting processing is sensitive 
not only to nonlinear scattering by microbubbles, but al- 
so to decorrelatton due to movement between pulses. 
Thirdly, since the bubbles are excited successively by 
interrogating pulses, they can be destroyed or in other 
ways altered between successive pulses. As a conse- 
quence, the resulting images, Doppler or other signals 
are not purely dependent on the nonlinear particularities 
of microbubbles, but also depend on these other possi- 
ble effects and artefacts. 

[0011] These techniques are sometimes designated 
" pulse-inversion imaging ", phase-inversion imaging ", 
"wideband harmonic imaging", "non-linear imaging us- 
ing phase cancellation etc. 

Summary of the invention 

[0012] The present invention pertains to a method of 
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coding the transmitted pressure-pulse waveform in such 
a way as to allow decoding the resulting echo wave- 
forms giving essentially zero contributions from linear 
scatterers or tissues which do not contain contrast agent 
microbubbles, and significant contributions for echoes 
originating from contrast agent microbubbles. 
[0013] The present invention relates to a method of 
nonlinear imaging by coding the transmitted waveforms 
within single transmit firings of the transducer or imaging 
probe. It is based on the observation that in the case of 
reflection by linear or essentially linear scatterers such 
as most tissues, individual rf-echo signals can be de- 
coded, or deconvolved, to produce essentially zero sig- 
nals, while in the case of reflection by more strongly non- 
linear scatterers such as contrast agent microbubbles, 
the same deconvolution algorithms produce significant- 
ly non-zero signals. 

[0014] In this way, a sensitive method is provided to 
significantly enhance the contrast between tissues and 
microbubbles, without paying a penalty of reducing the 
effective pulse-repetition frequency and frame rate. 
[0015] The invention also relates to a device for car- 
rying out the above-mentioned method, comprising 
means for constructing a double-pulse excitation wave- 
form, the two pulses composing this waveform having 
spectra with known relative frequency dependencies 
with respect to one another, differing in amplitude and 
phase in a known fashion, this waveform being intended 
to be applied to an ultrasound transducer to generate 
an ultrasonic beam, an ultrasound transducer array con- 
nected to said excitation means, comprising one or a 
plurality of transducer elements, a transmitter coupled 
to said transducer array for pulsing said transducer el- 
ements, receiving means coupled to said transducer for 
receiving said echo signals, means for deconvolving re- 
turned echo signal by a first appropriate decoding func- 
tion to obtain an rf-waveform, means for deconvolving 
returned echo signal by a second and different appro- 
priate decoding function to obtain an rf-waveform, 
means for realigning in time both deconvolved rf-wave- 
forms, means for normalizing amplitudes of said decon- 
volved rf-waveforms, and means for summing or other- 
wise combining both rf-waveforms to effectively cancel 
all echo components caused by linear scatterers. 
[0016] The invention further relates to the use of the 
device for detecting and imaging the nonlinear compo- 
nents of ultrasound echo signals returned from scatter- 
ers within the body of human patients or animals. 

Brief description of the drawing 

[0017] Figure 1 is a diagram of typical single-shot 
phase inversion pulse. 

[0018] Figure 2 is a diagram of amplitude spectra of 
single and double pulses. 

[001 9] Figure 3 is a diagram of transmit coding as del- 
ta function with respect to the first, "positive" pulse. 
[0020] Figure 4 is a diagram of transmit coding as del- 



ta function with respect to the second, "negative" pulse. 
[0021] Figure 5 is a diagram of of atypical rf-echo sig- 
nal. 

[0022] Figure 6 is a diagram of deconvolved rf-wave- 
5 form with respect to the "positive" pulse of figure 3. 
[0023] Figure 7 is a diagram of deconvolved rf-wave- 
form with respect to the "negative" pulse of figure 4. 
[0024] Figure 8 is a diagram of result of single-shot 
phase cancellation processing for a random collection 
10 of linear scatterers. 

[0025] Figure 9 is a diagram of deconvolved rf-wave- 
form according to pulse polarity 1 , with nonlinear scat- 
terers within a limited range. 

[0026] Figure 1 0 is a diagram of deconvolved rf-wave- 
15 form according to pulse polarity 2, with nonlinear scat- 
terers within a limited range. 

[0027] Figure 11 is a diagram result of single-shot 
phase cancellation processing for a random collection 
of nonlinear scatterers. 
20 [0028] Figure 12 is a block diagram which illustrates 
a contrast echography device with single-shot linear im- 
age cancellation imaging. 

Detailed description of the invention 

25 

[0029] In the disclosed implementation of this inven- 
tion, a pulse waveform such as illustrated in Fig 1 is gen- 
erated by a transmit ultrasound transducer, with a center 
frequency at 5 MHz. The horizontal time axis is labeled 

30 in microseconds, while the vertical axis is in arbitrary 
pressure or voltage units. This example illustrates a 
double-pulse, where the second instance is inverted in 
polarity, delayed by 0.5 u.S and reduced by 20% in am- 
plitude with respect to the first one. The double-pulse is 

35 constructed by summing two instances of the same 
model pulse, with the indicated changes in polarity, de- 
lay and amplitude; thus, the two pulses have spectra 
with identical relative frequency dependencies, differing 
only in amplitude and phase. The spectra magnitudes 

40 of the model waveform used and of the corresponding 
double-pulse waveform are given in Figure 2, illustrating 
that the double-pulse spectrum, when normalized in am- 
plitude to take into account its additional signal compo- 
nents, has as envelope the magnitude of the single- 

45 pulse spectrum, with a modulation given by the inverse 
of the time delay between the two pulses, i.e. 1/0.5 u.S 
= 2 MHz in this example. From linear-systems theory, 
the double-pulse waveform of Figure 1 can be consid- 
ered as a convolution of the model single-pulse with a 

so coding function composed of a pair of delta functions of 
amplitudes one and -0.8, separated by a time delay of 
0.5 u.S, as illustrated in Figure 3. 
[0030] Deconvolution of the received echo wave- 
forms from a collection of scatterers, such as those ob- 

55 tained in B-mode ultrasound imaging, can be applied in 
the time domain by convolution with an appropriate de- 
convolution kernel, or in the frequency domain by com- 
plex division by the spectrum of the coding of Figure 3. 
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One can thus obtain the rf-waveforms which would have 
been obtained from linear echoes to single-pulse exci- 
tation by the first pulse of the pulse-pair, by deconvolu- 
tion of the echoes from the pulse-pair by the appropriate 
decoding function. Similarly, one can also obtain the rf- 
waveforms which would have been obtained from linear 
echoes to single-pulse excitation by the second pulse 
of the pulse-pair, by deconvolution of the same echoes 
ffbm the pulse-pair by a different, appropriate, decoding 
function. 

[0031] This different decoding function can be ob- 
tained from the observation that the pulse pair can also 
be considered as a convolution of the second, negative 
phase, pulse component by a different coding of delta 
function pair, as given by Figure 4, where the first pulse 
is obtained by a delta function with a negative time shift 
and a 1/0.8, or 1 .25 amplitude relative to the delta func- 
tion of unity amplitude centered at the origin. In this way, 
and for any rf-echo signals produced by an ensemble of 
randomly distributed linear scatterers, it is possible, by 
deconvolution, to generate two rf-waveforms, equal to 
those that would have been obtained from the first and 
second pulses individually. Note that in this example, 
considering positive and negative pulses of different 
peak amplitudes allows computing the inverse of the 
coding spectra without any problems of discontinuities 
at the frequencies where zero amplitudes appear when 
the pulses are considered identical in amplitudes. 
[0032] Figure 5 illustrates the case of a typical rf-echo 
signal, for linear scatterers of random amplitudes and 
positions along the beam propagation path, in response 
to the "single-shot " double-pulse coded excitation. Two 
different deconvolved rf-waveforms can be derived from 
this signal : one obtained by deconvolution with respect 
to a first coding sequence, as given by Figure 3, and a 
second waveform obtained by deconvolution with re- 
spect to a second coding sequence, as given by Figure 
4. It is emphasized that both these resulting waveforms 
are derived from echoes originating from a single trans- 
mit, coded pulse waveform. The resulting rf-waveforms 
are illustrated in Figure 6 and Figure 7, respectively. The 
final step in the single-shot phase cancellation process- 
ing is to take the sum of both rf-waveforms, following a 
time shift to realign both rf-waveforms by an amount 
equal to the delay applied in the transmit double-pulse 
excitation waveform, and an amplitude normalization, 
again determined by that applied within the double- 
pulse sequence. Both these quantities can be accurate- 
ly known and programmed by the control circuits and 
software of modern echographic instruments. The result 
of that sum is an essentially zero time trace signal, as 
illustrated in Figure 8, indeed in agreement with what is 
predicted by linear systems theory for a collection of lin- 
ear scatterers. 

[0033] When the echo-signals are generated by non- 
linear scatterers as is the case in the presence of con- 
trast agent microbubbles following intravenous injection 
in the organism, the deconvolutions required to obtain 



traces such as those in Figure 6 and 7 no longer produce 
symmetrical results, and their sums produce non-zero 
output signals, revealing the presence of such non-lin- 
ear scatterers within the propagation path of the ultra- 

5 sound excitation beam. Localization of these waveforms 
is found from their times of arrival, just as in the case of 
standard pulse-echo processing. This situation is repre- 
sented in Figures 9 and 10, for deconvolutions accord- 
ing to pulse polarities 1 and 2, respectively, in the pres- 

10 ence of nonlinear scatterers at distances corresponding 
to times of arrivals between approximately 5.5 and 1 0.5 
|i8. 

[0034] Figure 1 1 illustrates the type of result that can 
be generated from the application of single-shot phase 

is cancellation processing along the lines described 
above. It can clearly be seen that the echoes arising 
from linear scatterers in regions outside 5.5 and 1 0.5 u,S 
are efficiently suppressed, while echoes arising from 
nonlinear scatterers within 5.5 and 1 0.5 \iS leave signif- 

20 icant non-zero signals available for further demodula- 
tion, processing and display as a representative of the 
presence of contrast agent in this region. 
[0035] In this way, a sensitive method is found to en- 
hance the contrast between regions not containing in 

25 significant amounts nonlinear scatterers, and regions 
containing contrast agent microbubbles with significant 
nonlinear scattering properties, without the drawbacks 
of other known methods, as outlined previously. It be- 
come evident from the disclosed method of the present 

30 invention that it relates to the field of imaging processing 
from ultrasound echoes reflected by scatterers, so that 
this method is neither a method for treatment of the hu- 
man or animal body, nor a diagnostic method practised 
on the human or animal body. Accordingly, this method 

35 is implemented by echographic imaging specialists and 
not by medical staff. 

[0036] A typical contrast echography device of the in- 
vention, with single-shot linear image cancellation for 
carrying out the above disclosed method is illustrated in 

40 figure 12. The device includes at least the following 
modules: central processor unit 1 , timing circuits 2, time 
gain control 3, transmit radiofrequency waveform shap- 
ing and coding 4, power transmit circuit 5, Tx/Rx-ele- 
ment multiplexer 6, ultrasound transducer array 7, re- 

45 ceive amplifiers with time gain function and analog to 
digital converter 8, digital beamforming processor 9, a 
first radiofrequency waveform deconvolver 1 0a, a sec- 
ond radiofrequency waveform deconvolver 10b, radiof- 
requency amplitude normalization and time shift module 

so 11a connected to deconvolver module 10a, radiofre- 
quency amplitude normalization and time shift module 
11b connected to deconvolver module 10b, digital sum- 
mation and one-line memory 1 2. The output of the mod- 
ule 12 is connected to a video scan converter 13, con- 

55 nected to a video monitor 14. 

[0037] In operation and under general control of a 
Central Processor Unit, the timing circuits 2 define a 
pulse repetition frequency,' required for constructing 
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echographic images (B-mode images, Doppler compo- 
nents in one- and two dimensions), based on sequential 
scanning of the region of the body to be imaged. For 
each successive coded-pulse excitation waveform, the 
timing circuits 2 also define the time-origin of a time de- 
pendent function used to provide variable amplification 
gain to echo signals originating from increasing imaging 
depths. This function is realized by the time gain control 
module 3, whose output can be a varying voltage, ap- 
plied to gain control of the receiving amplifier with ad- 
justable gain 8. The timing circuits 2 also allow to define 
transmit shaped/coded waveform, required for the ade- 
quate sequential excitation of the individual elements of 
the ultrasound transducer array 7 to provide beam fo- 
cusing and steering to be applied to the electrical exci- 
tation power transmit circuit 5. The timing circuits also 
provide the signals needed to bring predefined groups 
of array elements into connection with the power trans- 
mit circuit 5, by ways of the connections provided by the 
transmit-receive element multiplexer 6. The output sig- 
nals from this multiplexer are then routed to the receive 
amplifiers with time gain function and analog-to-digital 
conversion means 8, also controlled by the timing cir- 
cuits 2. The digital output of these circuits 8 is fed to the 
digital beamforming processor 9, the output of which is 
fed as a common input to the processing channels for 
implementing waveform decoding and linear image can- 
cellation on the returned echoes. 
[0038] The example of figure 2 implements linear can- 
cellation by feeding the echo signals through respective- 
ly first and second radiofrequency deconvolvers 1 0a, 
10b, each followed by a radiofrequency amplitude nor- 
malisation and time shift 1 1 a, 1 1 b in order to deconvolve 
the echo signals from the transmit shaped/coded exci- 
tation radiofrequency waveforms. The individual out- 
puts of the processing channels are then routed as input 
signals to the digital summation and one-line memory 
1 2 for taking the sum of both rf-waveforms to effectively 
cancel all echo components caused by linear scatterers 
and produce output signals then fed to the input of the 
video scan converter 1 3, setup to write the incoming da- 
ta, for each sequential pulse, in a pattern corresponding 
to the selected beam steering and positioning. 
[0039] Thus, by repeating the above sequence at the 
specified rate, each time modifying the beam steering 
and/or focusing to obtain echoes from successive posi- 
tions in the organs, tissue and blood, the scan converter 
output signal refreshes the image as displayed on the 
video monitor 14, in realtime, i.e. at: a rate between a 
few images per second to hundreds of images per sec- 
ond, sufficient for reproducing perception of movement 
by operator of the instrument. In the process described 
above, the regions of the echographic images corre- 
sponding to regions containing contrast agents appear 
with a contrast that is vastly enhanced, because tissues 
reflect minimal harmonic components. 
[0040] The benefits of the method and the system dis- 
closed may equally be exploited in systems in which the 



processing channels are part of the receiver of a one- 
or two-dimensional pulsed-Doppler ultrasound system, 
which may further incorporate a video output represent- 
ing a spectrum of velocity distribution and/or an audible 

5 signal output which is preferably a loudspeaker but it 
may also be any convenient sound reproducing device. 
Various useful options may be incorporated in the 
pulsed-Doppler ultrasound system such as a two-di- 
mensional map of velocity distribution which may further 

10 be colour coded, or it may incorporate a two-dimension- 
al map of echo amplitude or energy derived from Dop- 
pler echo components from moving targets, optionally 
at predetermined thresholds for velocities inferior or su- 
perior to a given value. 

is [0041] The deconvolving process may be implement- 
ed by spectral Fourier, Chirp-Z or wavelet transform 
analysis of returning echoes and also by applying the 
phase cancellation processing within a sliding time win- 
dow of the returning echoes. 

20 

Claims 

1. A method of detecting and imaging the nonlinear 
25 components of ultrasound echo signals returned 
from scatterers, characterized In that it comprises 
the steps of: 

a) constructing a double-pulse excitation wave- 
30 form, the two pulses composing this waveform 

having spectra with known relative frequency 
dependencies with respect to one another, dif- 
fering in amplitude and phase in a known fash- 
ion so that the double-pulse waveform can be 

35 considered as a convolution of any single of 

said pulses with a known coding function, the 
time delay separating the two pulses being 
shorter than the time interval needed for ultra- 
sound echo signals to be returned from the 

40 most distant said scatterers to be imaged, and 

longer than or equal to one half period of the 
fundamental central frequency of any single of 
said pulses, 

b) transmitting said double-pulse waveform 
45 within a single transmition of an ultrasound 

beam in the direction of said scatterers, 

c) receiving the echo signals in response to said 
transmitted beam, returned from said scatter- 
ers to form a received rf-waveform, 

so d) deconvolving said received rf-waveform by 

a first appropriate decoding function to obtain 
a new rf-waveform which would have been ob- 
tained from linear scatterers in response to sin- 
gle-pulse excitation by the first pulse of said 

55 double-pulse, 

e) deconvolving said received rf-waveform by 
a second and different appropriate decoding 
function to obtain another rf-waveform which 
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would have been obtained from linear scatter- 
ers in response to single-pulse excitation by the 
second pulse of said double-pulse, 

f) realigning in time both deconvolved rf-wave- 
forms obtained from the same returned echo s 
signal by a time delay as determined by the 
transmit-coding function, 

g) normalizing amplitudes of said deconvolved 
rf-waveforms, 

h) taking the sum of both rf-waveforms to effec- 1 o 
tively cancel all echo components caused by 
linear scatterers and produce output signals, 
and 

i) repeating the sequence b) to h) above in order 
to project ultrasound energy in a different direc- '5 
tion. 



2. A method of detecting and imaging the nonlinear 
components of an ultrasound signal as described in 
claim 1 , characterized In that the double-pulse ex- 
citation waveform is composed of two pulses having 
spectra with identical relative frequency dependen- 
cies and differing only in amplitude and phase in or- 
der to correspond to a convolution of a single of said 
pulses with a pair of delta functions of similar or 
slightly different amplitudes and separated by a 
known time delay. 

3. A method according to claim 1 or 2, wherein said 
returned echo signals are further demodulated and 
processed in known ways to allow their display on 
a video monitor as a two-dimensional map of echo- 
amplitudes. 

4. A method according to claims 1 -3, wherein a sum- 
mation algorithm is used for processing said output 
signals. 

5. A method according to claim 1 , wherein the ampli- 
tude of said output signals is coded by different vid- 
eo colours, to be superimposed on an otherwise 
conventional grey-scale video image obtained by 
the usual processing methods applied in B-mode 
imaging. 

6. A device for carrying out the method according to 
claim 1 , wherein it comprises: 

means for constructing a double-pulse excita- 
tion waveform, the two pulses composing this 
waveform having spectra with known relative 
frequency dependencies with respect to one 
another, differing in amplitude and phase in a 
known fashion, this waveform being intended 
to be applied to an ultrasound transducer to 
generate an ultrasonic beam, 
an ultrasound transducer array connected to 
said excitation means, comprising one or a plu- 



rality of transducer elements, 
a transmitter coupled to said transducer array 
for pulsing said transducer elements, 
receiving means coupled to said transducer for 
receiving said echo signals, 
means for deconvolving returned echo signal 
by a first appropriate decoding function to ob- 
tain an rf-waveform, 

means for deconvolving returned echo signal 
by a second and different appropriate decoding 
function to obtain an rf-waveform, 
means for realigning in time both deconvolved 
rf-waveforms, 

means for normalizing amplitudes of said de- 
convolved rf-waveforms, and 
means for summing both rf-waveforms to effec- 
tively cancel all echo components caused by 
linear scatterers. 

20 7. A device according to claim 6 comprising a video 
scan converter 

8. A device according to claim 6, comprising analog- 
to-digital converter circuits connected to said trans- 
25 ducer elements, whose purpose is to digitize the re- 
turned echo signals to allow processing of the re- 
turned echo signals by digital electronic processing 
circuits, and computing means. 

30 9, A device according to claim 6, wherein the returned 
echo signal processing means are part of a receiver 
of a pulsed-Doppler ultrasound system. 

1 0. A device according to claim 9, wherein said pulsed- 
35 Doppler ultrasound system incorporates an audible 

signal output by means of a loudspeaker. 

1 1 . A device according to claim 9, wherein said pulsed- 
Doppler ultrasound system incorporates a spectral 

40 video output representing a spectrum of velocity 
distribution. 

12. A device according to claim 9, wherein said pulsed- 
Doppler ultrasound system incorporates a two-di- 

45 mensional map of velocity distribution. 

13. A device according to claim 12, wherein said two- 
dimensional map of velocity distribution is colour 
coded. 

50 

14. A device claim 9, wherein said pulsed-Doppler ul- 
trasound system incorporates a two-dimensional 
map of echo-amplitude or energy derived from Dop- 
pler echo components from moving scatterers. 

55 

1 5. A device according to claim 9, wherein said pulsed- 
Doppler ultrasound system incorporates a two-di- 
mensional map of Doppler echo components from 
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scatterers moving with a velocity inferior to a pre- 
determined threshold. 

16. A device according to claim 9, wherein said pulsed- 
Doppler ultrasound system incorporates a two-di- $ 
mensional map of Doppler echo components from 
scatterers moving with a velocity superior to a pre- 
determined threshold. 

17. A device according to claim 6, wherein the required 10 
deconvolution processes of said returned echo sig- 
nals are obtained by spectral Fourier, Chirp-Z, 
wavelet transform analysis or direct time convolu- 
tion with a deconvolution function applied within a 
sliding time window on the returning echoes. 15 

18. A device according to claim 6, wherein said double- 
pulse waveform is constructed by a coding com- 
posed of a plurality of pulse counts higher than two, 
separated by known time delays and with known 20 
mutual spectral amplitude and phase relationships. 

19. Use of the device of one of claims 6-1 8 for detecting 
and imaging the nonlinear components of ultra- 
sound echo signals returned from scatterers within 25 
the body of human patients or animals. 
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